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ABSTRACT 

Context. In the absence of rotation and shear, and under the assumption of constant temperature or specific entropy, purely potential 
forcing by localized expansion waves is known to produce irrotational flows that have no vorticity. 

Aims. Here we study the production of vorticity under idealized conditions when there is rotation, shear, or baroclinicity, to address 
the problem of vorticity generation in the interstellar medium in a systematic fashion. 

Methods. We use three-dimensional periodic box numerical simulations to investigate the various effects in isolation. 
Results. We find that for slow rotation, vorticity production in an isothermal gas is small in the sense that the ratio of the root- 
mean-square values of vorticity and velocity is small compared with the wavenumber of the energy-carrying motions. For Coriolis 
numbers above a certain level, vorticity production saturates at a value where the aforementioned ratio becomes comparable with 
the wavenumber of the energy-carrying motions. Shear also raises the vorticity production, but no saturation is found. When the 
assumption of isothermality is dropped, there is significant vorticity production by the baroclinic term once the turbulence becomes 
supersonic. In galaxies, shear and rotation are estimated to be insufficient to produce significant amounts of vorticity, leaving therefore 
only the baroclinic term as the most favorable candidate. We also demonstrate vorticity production visually as a result of colliding 
shock fronts. 

Key words, magnetohydrodynamics (MHD) - turbulence - Galaxies: magnetic fields - ISM: bubbles 



o 



^ ! 1- Introduction 



CZ5 
C3 



(N 
> 

oo 

(N 

in 

oo 
o 

O 



X 



Turbulence in the interstellar medium (ISM) is believed to be 
driven by supernova explosions. Such events inject sufficient 
amounts of energy to sustain turbulence with rms velocities of 
~ lO km/s and correlation lengths of up to 100 pc dBeck et all 
Simulations of such events can be computationally quite 
demanding, because the bulk motions tend to be supersonic and 
the flows involve strong shocks in the vicinity of individual ex- 
plosi on sites, as was seen early on in two-dimensional simula- 
tions (Ros en & BregmanLll995h . Nevertheless, such simulations 
are able to reproduce a number of physical phenomena such 
as the observed vo l ume fractions of ho t, warm, and cold gas 
dRosen et all Il996t iKorpi et all I l999al) . the statistics of pres- 
sure fluctuations ( Mac Low et all 1 2005). the effects of the mag- 
netic fi eld dde Avillez & Breitschwerdt. 2005), an d even dynamo 
action iGressel et aL, 2008t lGissinger et al., 2009; Hanasz ^t all 
2009). These simulations tend to show the development of sig- 
nificant amounts of vorticity, which is at first glance surprising. 
Indeed, each supernova drives the gas radially outward and can 
roughly be described by radial expansion waves. In such a de- 
scription, turbulence is forced by the gradient of a potential that 
consists of a time-dependent spherical blob at random locations. 
Obviously, such a forcing is irrotational, so no vorticity is pro- 
duced. 



Earlier work of Mee & Brandenburg (2006) showed that un- 
der isothermal conditions only the viscous force can produce 
vorticity and that this becomes negligible in the limit of large 
Reynolds numbers or small viscosity. In principle, vorticity can 
also be amplified akin to the dynamo effect by the V x (u x ui) 
term, which is analogous to the induction term in dynamo the- 
ory, where u) plays the role of the magnetic field. However, nei- 



ther this nor the viscosity effect were found to operate - even 
at numerical resolutions of up to 512 3 meshpoints. This dis - 
agreed with subsequent simulations by Federra th et all (|2010), 
who solved the isothermal inviscid Euler equations with irrota- 
tional forcing using the Flash Code. They found significant 
vorticity generation in proximity to shocks where some kind of 
effective numerical viscosity must have acted. 

Given that under isothermal conditions, only viscosity can 
lead to vorticity production, one must ask whether numerical vis- 
cosity or effective viscosity needed to stabilize numerical codes 
might have contributed to the production of vorticity in some of 
the earlier works. Indeed, it is possible that the directional opera- 
tor splitting used in the Flash Code may have b een responsible 
for spu rious vorticity generation in the work of Federrath et all 
(2010); (R. Rosner, private communication). On the other hand, 
when cooling and heating functions are included to perform 
more realistic simulations of the ISM, vorticity could be pro- 
duced by the baroclinic term. Furthermore, even in the isother- 
mal case, in which the baroclinic term vanishes, vorticity could 
be produced if there is rotation and/or shear. 

The baroclinic term results from taking the curl of the pres- 
sure gradient term and is proportional to the cross product of the 
gradients of pressure and density. This term can play an impor- 
tant role when the assumptions of isothermality or adiabaticity 
are relaxed. Indeed, the baroclinic term can also be written as the 
cross product of the gradients of entropy and temperature. This 
formulation highlights the need for non-ideal effects, because 
in the absence of any other heating or cooling mechanisms, the 
entropy is just driven by viscosity. Again, it is not obvious that 
in the absence of additional heating and cooling much vorticity 
can be produced. On the other hand, it is clear that viscous heat- 
ing must be significant even in the limit of vanishing viscosity, 
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because the velocity gradients can be very large, especially in 
shocks. Of course, the assumption about additional heating and 
cooling is not realistic for the interstellar medium and will need 
to be relaxed. Finally, there are the effects of rotation and shear, 
that can contribute to the production of vorticity even in the ab- 
sence of baroclinicity. 

The goal of this paper is to study the relative importance of 
the individual effects that contribute to vorticity production. It is 
then advantageous to restrict oneself to simplifying conditions 
that allow one to identify the governing effects. An important 
simplification is the restriction to weakly supersonic conditions 
so that shocks and other sharp structures can still be resolved 
with just a uniform and constant viscosity. We also neglect the 
effects of stratification which can only indirectly contribute to 
vorticity production. In fact, a constant gravitational accelera- 
tion drops out when taking the curl. Only in the non-isothermal 
and non-isentropic case can gravity contribute to vorticity pro- 
duction by enhancing the effect of the baroclinic term. We begin 
with a preliminary discussion and a qualitative analysis of the 
important terms in the vorticity equation. 

2. Preliminary considerations 

We recall that in the absence of baroclinicity, rotation, and shear, 
the curl of the evolution equation of the velocity is given by (see, 
e.g JMee & Brandenburg! 120061) 

du) 

— = V x (it x u> - i>V x u) + vV x G, (1) 

where v is the kinematic viscosity (assumed constant) and Gi — 
2Sij Vj In p is a part of the viscous force that has non-vanishing 
curl even when the flow is purely irrotational. Here, 

Sij = \{u it j + Uj ti ) - \5ijU ktk (2) 

is the traceless rate of strain matrix, and commas denote par- 
tial differentiation. The G term breaks the formal analogy with 
the induction equation. It is convenient to express the resulting 
rms vorticity in terms of the typical wavenumber fc w of vortical 
structures which we define as 

^cj ^rms/^rms- (3) 

We monitor the ratio fc ^/fcf, where fcf is the adopt ed nominal 
forcing wavenumber. In lMee &~Br andenburg (2006), the result- 
ing vorticity, expressed in terms of the ratio k u /k{, was found to 
be zero within error bars. This result is compatible with the idea 
that the v V x G term in Equation (Q3 is insignificant for vorticity 
production. By contrast, in vortical turbulence and at moderate 
values of t he Reynolds numb er, k^/ki is found to be of the order 
of unity ( Brandenb urg[l200"ll) . although one should expect a mild 
increase proportional to the square root of the Reynolds number 
as this number increases. 

2.1. Rotation 

Rotation leads to the addition of the Coriolis force, 2f2 x u, 
in the evolution equation for the velocity. Taking the curl, we 
obtain the vorticity equation (Q} with two additional terms, both 
proportional to so we have 

— = ...- 2nv± ■ u ± + 2n ■ Vm_l, (4) 

ot 

where the dots denote the other terms in Equation (Q]i that we dis- 
cussed already. In order to estimate the production of vorticity, 



one could derive an evolution equation for the enstrophy density, 
t;u> 2 , by multiplying the right-hand side of Equations (Q} and (0]i 
by ui, and use a closure assumption for the resulting triple corre- 
lations. However, it is then difficult to obtain a useful prediction 
for w rms , because the right-hand side of such an equation would 
necessarily be proportional to u> and would therefore vanish, un- 
less u) lma was different from zero to begin with. Instead, we esti- 
mate w rms by computing the rms value of du> jdt and replacing 
it by w rm s/Tsi, where tq is a typical time scale of the problem. 
This leads to 

co rms « 2Qtq ((V± • u ± ) 2 + (V||« ± ) 2 ) 1/2 , (5) 

where Vj_ and V|| denote derivatives in the directions perpen- 
dicular and parallel to the rotation axis and u± is the velocity 
vector perpendicular to the rotation axis. Using Cartesian coor- 
dinates where CI points in the z direction, we have 

1/2 

LJ Ims « 2Qtq ((u x , x + u V: y) 2 + u x z + u VtZ ) . (6) 

We expect tq to be comparable to the turnover time, r = 
(ttrmsfcf) -1 . We expect the rms values of the velocity deriva- 
tive term in Equation ((6} to be comparable to the rms velocity 
and some inverse length scale. Typically, one would expect it to 
be proportional to u lms k{, although, again, there can be an ad- 
ditional dependence on the square root of the Reynolds number. 
However, for fixed Reynolds number, and not too rapid rotation, 
we expect w rms to increase linearly with the Coriolis number, 
i.e., 

Co = 2Qt, where r = (ur ms k{)~ . (7) 

Thus, we expect fc w /fcf = St^Co, where we have defined an 
effective rotational Strouhal number, 

St n = T£ ff u rms fc f . (8) 

We regard this as a fit parameter that will emerge as a result 
of the simulations. We have here introduced the quantity r n s , 
where Tq H /tq is given by the ratio of the velocity gradient terms 
divided by u rms fcf . However, for larger values of Co there may 
be a departure from a linear dependence between k^/ki and Co. 
(We note that, apart from a possible 4ir factor, the Coriolis num- 
ber is just the inverse Rossby number.) One aim of this paper 
is therefore to verify this dependence from simulations and to 
determine empirically the value of tq. 

2.2. Shear 

In the presence of linear shear with u s = (0, Sx, 0), the evolu- 
tion equation for the departure from the mean shear attains addi- 
tional terms, — w s V ■ u — u ■ Vu s . This implies a dependence 
of w rms on S, analogous to the 57 dependence discussed above. 
In components form, this means that 

1/2 

Wrms ~ St S ((U x<x + Uy t y) 2 + U x z + u\ y + 0(xu")) , (9) 

which is quite similar to Equation (0, except that in the penulti- 
mate term in angular brackets the indices are now interchanged, 
i.e. we now have u z , y instead of u y . z . In analogy to tq, we de- 
fine ts as a typical time scale of the problem and we expect 
it to be again related to the turnover time r. The 0(xu") de- 
notes the presence of additional terms that are proportional to 
x and to second derivatives of u. However, when adopting the 
shearing box approximation with shearing-periodic boundarie s 
dGoldreich & Lvnden-BellL 119651: IWisdom & TremaineL [l988), 
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each point in the xy plane is statistically equivalent. We would 
therefore not expect there to be a systematic x dependence, 
which corresponds to the assumption of Galilean invariance that 
is sometimes used in t he study of turbulent tran s port c oefficients 
in linear shear flows (Sridhar & Subramanian, 2009). We will 
postpone the possibility of additional terms until later. Since we 
expect ts to be comparable to r = (u rms fcf the rms vorticity 
should be proportional to the shear parameter, 



Sh = St = S/u rms k f , 



(10) 



although for large values of |Sh| we may expect departures from 
a linear dependence. Determining this dependence is another 
aim of this paper. Again, a linear dependence is characterized 
by the values of t$ and r£ ff , where, in analogy with the previ- 
ous case with rotation, the ratio t$ s /ts is given by the deriva- 
tive term in Equation (|9), normalized by u lms k{. A convenient 
non-dimensional measure of the value of t s s is what we call the 
shear Strouhal number, 

St S = rfurmakt, (H) 

which can be determined provided there is a range in Sh over 
which u> lms increases linearly with Sh. 

The study of vorticity production by rotation and shear is 
quite independent of thermodynamics and can in principle be 
studied even in the incompressible case. However, in the present 
paper we study this effect in the weakly compressible case of 
low Mach numbers and under the assumption of an isothermal 
equation of state, where the baroclinic term vanishes. 

2.3. Baroclinicity 

As mentioned in the introduction, the baroclinic term, propor- 
tional to Vp x Vp, emerges when taking the curl of the pressure 
gradient term, p -1 Vp. This term can also be written as 



p- 1 Vp = Vh- TVs, 



(12) 



where h and s are specific enthalpy and specific entropy, respec- 
tively, and T is the temperature. Thus, we have 



— = ... + VTx Vs. 

ot 



(13) 



In order to study the effect of the baroclinic term, it is useful to 
look at the dependence of the mean angle 9 between the gradi- 
ents of s and T, defined via 



sin 2 d = ((VT x Vs) 2 )/((VT) 2 ){(Vs) 2 ). 



(14) 



An important aspect is then to study first the dependence of the 
rms values of the gradients of s and T. We can do this by looking 
at a one-dimensional model where, of course, 6 — 0. 

Next, we need to determine 9 from three-dimensional simu- 
lations. The hope is then that we can express baroclinic vorticity 
production in the form 



k u /k{ = Stbaro(VT) rms (Vs) rms sin6'/u 2 ms fc f 2 



On dimensional grounds we expect the product of ( VT) rms and 
(Vs) lms to be of the order of u 2 ms k 2 , and so a possible ansatz 
would be 



k u /kf = Stbf ro sin0 



(16) 



where we have subsumed the scalings of (VT) rms and ( Vs) r 
in that of an effective baroclinic Strouhal number St£f ro . 



An important issue is the fact that viscous heating leads 
to a continuous increases of the temperature. As a result, the 
sound speed changes and it becomes then impossible to study 
the behavior of the system in a steady state. In order to avoid 
this inconvenience, we add a volume cooling term that is non- 
vanishing when the local sound speed c s is different from a given 
target value, c s o- Thus, in the presence of finite thermal diffusiv- 
ity x, and with a cooling term governed by a cooling time r coo i, 
our entropy equation takes the form 



T^£ = 2 ,S 2 
Dt 



1 



p" 1 V • (c p pxVT) (cj - c 2 s0 ), (17) 



7"cool 



where c s is the adiabatic sound speed. We assume a perfect gas 
so that c 2 = (7 — l)c p T, where 7 = c p /c v = 5/3 for a 
monatomic gas, and c p and c v are the specific heats at constant 
pressure and constant volume, respectively. The value of r coo i 
can have an influence on the results, so we need to consider dif- 
ferent values. We express r coo i in terms of c s o and fcf, and define 
the nondimensional quantity St coo i = r coo ic s ofcf . 

3. The model 

In this paper we solve the continuity equation for the density p, 
Dlnp 



Dt 



= -V • u, 



together with the momentum equation for the velocity it, 
Dm 

— = -p -1 Vp -20xu- Su x y + V<f> + F visc , 



(18) 



(19) 



where D/Di = d/dt + (u + u s ) ■ V is the advection operator 
with respect to the sum of turbulent flow u and laminar shear 
flow u s , p is the pressure, <j> is the forcing potential, and 



p- J V • (2vpS) 



(20) 



is the viscous force, where S was defined in Equation (f2]). The 
forcing potential is given by 



(x,t) = 0o Ncxp{-[x - Xi{t)] 2 /R 2 } 



(21) 



where x = (x, y, z) is the position vector, Xf(t) is the random 
forcing position that changes abruptly after a time interval At, R 
is the radius of the Gaussian, and N is a non-dimensional factor 
proportional to At' 1 / 2 . This ensures that the amplitude of the 
correlation function of <fi is independent of At. Thus, we choose 
N = y/R/c s oAt. Since N is non-dimensional, the prefactor 0o 
has the same dimension as <fi, which is that of velocity squared. 
We consider two forms for the time dependence of Xf . First, we 
take Xf such that the forcing is <5-correlated in time. In that case, 
At is equal to the length of the time step St. Alternatively, we 
choose a finite forcing time Stf OICC that defines the interval dur- 
ing which X{ remains constant, after which the forcing changes 
again abruptly. Thus, 



(15) At = max (St, St loicc ) 



(22) 



is equal to St in the ^-correlated case or equal to 8tf OICC in the 

case of finite corr elation time. 

The work of Mee & Brandenbur 3 d2006h showed that the 
peak of the energy spectrum depends on the radius R of the 
Gaussian. Indeed, the Fourier transform of exp(— r 2 /R 2 ) is also 
a Gaussian with cxp(— k 2 /kf), where 



2/R. 



(23) 
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In the following we use this as our definition of fcf and check a 
posteriori that this is close to the position of the peak of the ener- 
gy spectrum. In the following, we characterize our simulations 
in terms of the ratio fcf/fci, and consider values between 2 and 
10. 

We use the PENCIL CodeQ which is a non-conservative, 
high-order, finite-difference code (sixth order in space and third 
order in time) for solving the compressible hydrodynamic and 
hydromagnetic equations. We adopt non-dimensional variables 
by measuring speed in units of a reference sound speed, c s q, and 
length in units of 1/fci, where fci is the smallest wavenumber in 
the periodic domain. This implies that the nondimensional size 
of the domain is (27r) 3 . 

In order to study the effects of rotation and shear, we ignore 
entropy effects and restrict ourselves to an isothermal equation 
of state with constant sound speed c s . This means that Vp 
reduces to cjVlnp = V/i, which has vanishing curl. Here, 
h = c 2 In p is the relevant enthalpy in the isothermal case. On 
the other hand, in order to study the effects of baroclinicity, we 
do need to allow the entropy to vary, so we also need to solve 
Equation ( fTTI ). and study the dependence of k u /k{ on the Mach 
number, 



Ma = u lms /c s 



(24) 



In order to characterize the degree of turbulence, we define the 
Reynolds number based on the energy-carrying scale, corre- 
sponding to the typical wavenumber where the spectrum peaks, 
i.e. 



Re = u Yms /uk f . 



(25) 



For vortical turbulence, this definition is known to be a good 
measure of the ratio of the r esulting turbulent v iscosity divided 
by the molecular diffusivity dYousef et aU l2~003). The two num- 
bers, Ma and Re, can be varied by changing v and/or the strength 
of the forcing. In all cases we use \ — v - Another input param- 
eter is the forcing Strouhal number 



Stfo 



(26) 



which is zero for (5-correlated forcing and equal to about 0.3 in 
ca ses with finite correlation tim e. These are also the values used 
bv lMee & Bran denburg (200(|. 

In the following we also consider kinetic energy and enstro- 
phy spectra, Ex(k) and E u ,( k), respectively. They are normal- 

— m> 



ized such that (lLesieurill99C 



E K (k)dk = ±(m 



E u {k)dk = \(u 2 



(27) 



where I 



u 2 ) and \{u 2 ) are kinetic energy and enstrophy, re- 
spectively. For comparison we also consider spectra of enthalpy, 
Eh(k), which are normalized such that J E^(k) dk = \(h 2 ). 

Throughout this paper we assume periodic boundary condi- 
tions, except that in the presence of shear we employ shearing- 
periodic boundary conditions where the x direction is periodic 
with respect to positions in y that shift with time, i.e. 



f(-hL x ,y, z,t) = f(±L x ,y + L x St, z,t), 



(28) 



where / represents any one of our four dependent vari- 
ables (u,p). This bound a ry co ndition was first proposed by 
iGoldreich & Lvnden-Belll ( 1 1965b and has been routinely used 
in local simulations of accretion disk turbulence (Hawlev et al., 
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Fig. 1. Time-averaged kinetic energy and enthalpy spectra for 
two values of the Coriolis number for Re = 25 and Stf orce = 0.4. 
The two straight lines give the slopes —2 and —3, respectively. 
In both cases we have fcf /fci = 4. 
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Fig. 2. Dependence of fc w /fcf on Co for three groups of runs: 
group 1 with Re = 15, fcf/fci = 10, Stf orce = 0.3; group 2 with 
Re between 25 and 130, fcf/fci = 4, Stf orcc = 0.4; and group 3 
with Re = 30 and 150, fcf/fci = 2, St forcc = 0. 



1 19951). Not e, however, t hat re cent work of iRegev & Um urhan 



(2008) and iBodo et ail (120081) called attention to the possibil- 
ity of problems with the shearing sheet approximation when the 
size of the perturbations is large. In somewhat weaker form, this 
problem also applies to a non-shearing periodic box. Indeed, we 
shall kept this in mind when interpreting some of the results pre- 
sented below. 



4. Results 

We begin by studying the effect of rotation. In Figure Q] 
we plot time-averaged kinetic energy and enthalpy spectra, 
E}<i{k) and Eh(k), respectively. Note that rotation has a ten- 
dency to move the peak of E-&(k) to the left of the nom- 
inal value of fcf. However, at the Reynolds number of 25 
shown here, there is no inertial range, but in all cases, the en- 
ergy spectra show a clear viscous dissipation range, suggest- 
ing that these runs are sufficiently well resolved. At somewhat 
larger Reyn olds number or smal l er for cing wavenumber, ear- 
lier work of iMee & Brandenburg! (120061) began to show a short 
k~ 2 subrange. Such a slope is predicted for shock turbulence 
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Fig. 3. Dependence of fc w /fcf on Re for Co w 1, fef/fei = 2, and 
Stforce = (^-correlated forcing). 
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Fig. 4. Time-averaged enstrophy spectra, E u (fc), compared with 
k 2 EK{k), for Re = 25, Stf orce = 0.4, and three values of the 
Coriolis number. The curves of k 2 EK(k) are close together and 
overlap for Co = 0.01 (dotted) and 0.15 (dashed), so it becomes 
a single dash-dotted line. The fc -3 slope is shown for compari- 
son. In all three cases we have fcf/fci = 4. 



dKado mtsev & Petviashvijjl Il973h . and it has also been seen in 
the ir rotational component of transonic turbulence dPorter et al.L 
[l998h . 

In Figure [2] we plot the dependence of fc w /fcf on Co for 
three groups of runs: group 1 with Re = 15, fcf/fci = 10, 
Stforco = 0.3; group 2 with Re between 25 and 130, fcf/fci = 4, 
Stforco = 0.4; and group 3 with Re = 30 and 150, fcf/fci = 2, 
Stf orco = 0. In all cases we use 128 3 mesh points, average the 
results over at least 200 turnover times and, in some cases, even 
several thousand turnover times. It turns out that for Stf orC c / 
a linear relationship between fc^/fcf and Co is a good approxi- 
mation for Co <; 10, where k u /kf » 0.03 Co, i.e. St n = 0.03. 
Furthermore, we see from Table Q] that the normalized veloc- 
ity derivative terms are all about 0.5, so the root of the sum 
of their squares is slightly larger than unity, corresponding to 
Tq r/rn ~ 1.3. For Co > 10 the value of fc w /fcf seems to satu- 
rate at about unity. 

A similar result is also found for Stf orce = 0, except that 
there remains a spurious contamination of vorticity even for 
small values of Co, a limit in which we expect to observe no 
vorticity production. By varying the value of Re, while keep- 
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Fig. 5. Dependence of fe^/fcf on Sh for Re w 40 and fcf/fci = 2 
and 5-correlated forcing (Stf orcc = 0). Different resolutions are 
shown to give similar results. At small values of |Sh|, compar- 
isons with Stf orC o = 0.3 (keeping fcf/fci = 2) or fcf/fci = 10 
(keeping Stf orC c = 0) are also shown. 



ing Co w 1 fixed, we see that k^/kf asymptotes to zero for 
sufficiently small values of Re; see Figure [3] This suggests that 
there can easily be spurious vorticity generation, possibly due 
to marginally sufficient resolution. The possibility of spurious 
vorticity is indeed verified by Figure [4] where we compare en- 
strophy spectra, E u (k), with fc 2 _Ex(fc)- Note that for large val- 
ues of Co, the enstrophy spectrum decays like fc~ 3 . However, for 
smaller values of Co the level of enstrophy at the mesh scale re- 
mains approximately unchanged and is thus responsible for the 
spurious vorticity found above for small values of Co and not 
too small values of Re. Nevertheless for larger values of Co, the 
production of vorticity is an obvious effect of rotation in an oth- 
erwise potential velocity field, and it is most pronounced at large 
length scales, as can also be seen in Figure [4] 

Next, we study the dependence of the ratio k^/ki on shear; 
see Figure|5] We use a resolution of 64 3 or 128 3 mesh points, av- 
erage the results over at least 200 turnover times and, in cases of 
lower resolution, over several thousand turnover times. It turns 
out that in the presence of shear, some level of helicity pro- 
duction can never be avoided - even in the limit of small Sh. 
Again, this appears spurious and demonstrates the general sen- 
sitivity of vorticity generation on resolution effects. An addi- 
tional problem is of course the finite size of the shearing box 



Table 1. Root-mean-squared values of components of the veloc- 
ity derivative tensor, normalized by M rm sfcf , as well as the three 
diagonal components of the (uiUj) tensor for 4 values of Co. 



Co 


0.11 


0.35 


0.99 


2.80 


(V±u±) rms /u Tms kf 


1.26 


1.20 


1.21 


1.04 


l^x^x /^rms^f 


0.76 


0.74 


0.75 


0.70 


,y /^rms^f 


0.79 


0.74 


0.75 


0.70 




0.49 


0.47 


0.48 


0.63 


7/ rms hi hr 
U'yjZ 1 "Tins "'i 


0.49 


0.47 


0.48 


0.63 


^>z,y /Wnns^f 
i^x) 1 ^rms 


0.49 


0.47 


0.46 


0.58 


0.32 


0.33 


0.35 


0.42 


(^|//^rms 
{^2 ) /^rms 


0.34 


0.34 


0.35 


0.41 


0.35 


0.34 


0.29 


0.18 
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Fig. 6. Dependence of the rms values of temperature and entropy 
on </>o for v/c s R = 1 and St coo i = 0.2 (top panel), 0.6 (middle 
panel), and 2 (bottom panel). 



dRegev & Umurhanl 120081: iBodo et al.Ll2008l) . which may be re- 
sponsible for spurious vorticity generation. On the other hand, 
there is vorticity generation even for large scale-separation ra- 
tios, kf/ki = 10; see the dash-dotted line in Figure [5] This 
suggests the possibility of a more general problem that would 
not go away even in the limit of small eddies and small val- 
ues of |Sh|. Nevertheless, there is a clear rise of k u /k{ when 
|Sh| > 0.1, which is in agreement with our expectations out- 
lined in Section l2~2l However, the slope in this relationship is 
rather steep, Sts w 6. The velocity derivative terms are only 
slightly larger than in the case with rotation, corresponding to 
T S l T s ~ 1-5; see a ls° Table [2] Tentatively, this suggests that 
for comparable values of Co and Sh, rg 3> tq. On the other 
hand, given that even for small values of Sh there is spurious 
vorticity generation, we cannot be certain that the results are re- 
liable for larger ones either. The case with shear must therefore 
remain somewhat inconclusive. 



Table 2. Similar to Table [TJ but for the case with shear. 



Sh 


-0.01 


-0.03 


-0.06 


-0.12 


-0.26 


(V_Lttx)rms/Urmsfcf 


1.35 


1.36 


1.36 


1.37 


0.87 


^x,x /^rms^f 


0.87 


0.87 


0.90 


0.97 


0.75 




0.81 


0.81 


0.76 


0.66 


0.47 


^x,z /Wrms^f 


0.51 


0.50 


0.51 


0.53 


0.72 


^y,z /^rms^f 


0.46 


0.47 


0.47 


0.46 


0.74 


(^}/^rms 


0.46 


0.48 


0.46 


0.43 


0.57 


0.37 


0.38 


0.37 


0.36 


0.25 


(t/^) /l£ rms 


0.31 


0.31 


0.32 


0.34 


0.56 




0.32 


0.31 


0.31 


0.31 


0.25 
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Fig. 7. Dependence of Mach number and the rms value of en- 
tropy on <fio for v/c s R = 1 and St coo i = 0.2, 0.6, and 2 (solid, 
dotted, and dashed line types, respectively). 



Finally, we consider the possibility of vorticity generation by 
the baroclinic term. In a preparatory step we study first the de- 
pendence of the product (VT) rms (Vs) rms on both Ma and Re in 
a one-dimensional model. In all cases we vary the strength of the 
forcing amplitude in the range 1 < 4>o/c^ Q < 500 for different 
values of viscosity and cooling time. As we increase the value of 
0o, the Reynolds number increases for a given value of the vis- 
cosity. For small values of <fio, the Mach number also increases 
linearly, where the ratio of Ma/Re increases with increasing vis- 
cosity. However, for larger values of </>o there is saturation and 
Ma no longer increase with <f>Q. 

Furthermore, in the range where Ma still increases linearly 
with <fio, the rms value of the entropy gradient increases, but it 
also saturates when Ma saturates. The rms value of the temper- 
ature gradient, however, decreases with increasing values of </>o, 
but this seems to be a special property of the one-dimensional 
model that is not borne out by the three-dimensional simulations 
where it stays approximately constant. 

Remarkably, the results are fairly independent of the cooling 
time, except that the breakpoint where (Vs) rms saturates occurs 
for smaller values of </>o as we increase the cooling time; see 
Figure [6] This break point is also related to the point where the 
Mach number saturates, as can be seen from Figure|7] However, 
for longer cooling times there can be longer transients, mak- 
ing it difficult to obtain good averages. Therefore we focus, in 
the rest of this paper, on the case of shorter cooling times using 
Stcooi = 0.2. Another remarkable result is that the normalized 
value of (VT) rms ( Vs) rms is always of the order of about 10 -3 , 
independent of resolution, cooling time, or the value of the vis- 
cosity. 

Most of the basic features of the one-dimensional model 
are also reproduced by two- and three-dimensional calculations. 
Two-dimensional simulations have the advantage of being easily 
visualized and are therefore best suited for illustrating vorticity 
production by the baroclinic term. In Figure [8] we demonstrate 
that vorticity production is associated with the interaction be- 
tween the fronts of different expansion waves. In this example 
we chose 8tf mce c s o/R = 0.1, so the first expansion wave is 
launched at t = and the second one at tc a o/R = 0.1. The 
top row of Figure[8]shows that at tc s o/R = 0.09, i.e. just before 
launching the second expansion wave, the baroclinic term and 
the vorticity are still just at the noise level. At that time the most 
pronounced feature is the discontinuity between the Gaussian 
expansion waves in the periodic domain. This leads to negligi- 
bly weak vorticity, and no baroclinic term. At tc a o/R = 0.11, 
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Fig. 8. Images of T, s, (VT x V.s) z , and normalized vertical vorticity for a two-dimensional run with 5tf mce c s o/R = 0.1 at an 
instant shortly before the second expansion wave is launched (top row), and shortly after the second expansion wave is launched 
(second and third row). Note the vorticity production from the baroclinic term in the second and third row, while in the top row, 
(VT x Vs) z and u z are just at the noise level of the calculation. Even under our weakly supersonic conditions shock surfaces are 
well localized and the zones of maximum production of vorticity appear to be those in which the fronts encounter each other. Here 
we have used <fio / c 2 — 100, v = x — O.lcgo-R, with 512 2 mesh points. Only the inner part of the domain is shown. 



the effect of the second expansion wave becomes noticeable in 
visualizations of both (VT x V.s) 2 and lu z , while our visual- 
izations of T and s barely show the second expansion wave. At 
tc s o/R = 0.14, the first expansion wave is clearly no longer cir- 
cular, which is obviously associated with the second expansion 
wave that is now quite pronounced in our visualizations of both 
Tands. 

In order to have a more accurate quantitative determina- 
tion of vorticity production, we now consider three-dimensional 
models. In Figure [9] we show the dependence of various quan- 
tities on (f>o for St coo i = 0.2 and v/c s R = 1. In all cases we 
use 128 3 mesh points and average the results over between 20 
and 70 turnover times. Note that here Re sa 0.05 4>q/c 2 q . Given 
that Re depends inverse proportionally on v/Rc s $, we can also 
write Re w 0.05 4)qR/c s qv. The Mach number saturates at about 
Ma = 3, and the rms value of the entropy gradient increases up 
until this point. Given that the rms value of the temperature gra- 
dient also stays approximately constant, we find a weak increase 
of (VT) rms (Vs) rms . The value of (VT x Vs) rms is always 
found to be a certain fraction below this value, resulting in a 
typical baroclinic angle of about 45 degrees; see the third panel 



of Figure [9] Finally, the amount of vorticity production in terms 
of k u /k{ is about 0.3 for (j>o/c 2 <; 20. For smaller values, on 
the other hand, there is an approximately linear increase with 
Wfcf~ 0.014 o /c 2 o . 

The possibility of spurious vorticity is easily eliminated in 
this case by looking at enstrophy spectra; see Figure [10] where 
we compare E LJ (k) with k 2 E-^(k). All spectra fall off rapidly 
with increasing k. Thus, even though the initial vorticity genera- 
tion occurred evidently at the smallest available scales, once the 
flow becomes fully developed, most of the enstrophy resides at 
scales equal to or larger than the driving scale. Furthermore, the 
spectra of E u (k) and k 2 EK(k) are close together, suggesting 
that the vorticity is close to its maximal value. 



5. Applications 

The level of vorticity that is produced in the usual case of 
solen oidal forcing of the turbulence is such that k u /ki ps 1 (see, 
e.g., [Brandenburg] 1200 ll) . For turbulence whose driving force 
has finite correlation time (Stf orcc = 0.3, for example), and small 
values of Re, we have k u /k{ = 0(1) when Co ^ 10. However, 
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Fig. 9. Dependence of Ma, Re, the rms values of Vs and 
VT, the angle 8 between them, as well as k^/kf, on 4>o/c 2 for 
v/c s R = 1. 



for larger values of Re, the turbulence becomes vortical already 
for smaller values of Co. Comparing with the galaxy, we have 
£1 k 10 5 s , Urms = lOkm/s, and an estimated correlation 
length of about 70 pc, so k f = 3 x 10~ 20 cm, so Co = 0.07, 
which is rather small. Thus, rotation may not be able to pro- 
duce sufficient levels of vorticity. Given that in galaxies with flat 
rotation curves, S « — fi, shear should not be very efficient ei- 
ther. However, the Mach numbers are undoubtedly larger than 
unity in the interstellar medium, so this should lead to values of 
ku/kf w 0.3, which is the saturation value found in Figure [9] 
Given that one of the reasons for studying the production of vor- 
ticity is the question of dynamo action, we should point out that 
such values of k u /k{ are large enough for the small-scale dy- 
namo. Large-scale dynamo action should be possible in galax- 
ies as well, because of their large length scales, but it suffers 
from the well-known problem of a small growth rate. It then 
remains difficul t to explain large -scale magnetic fields in very 
young galaxies (Beck et aTL ll996h . 
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Fig. 10. Time-averaged enstrophy spectra, E a) (k) (thick lines), 
compared with k 2 E\^(k) (thin lines below the correspond- 
ing thick lines), for the three-dimensional baroclinic case with 
O / C 2 O = io (dashed), 100 (dotted), and 500 (solid lines). The 
k~ 2 slope is shown for comparison. In all three cases we have 
k { /ki = 4. 



The question of vorticity generation is also important in stud- 
ies of the very early Universe, where phase transition bubbles 
are believed to b e generated in connection with the electroweak 
phase transition (Kaianti e & Kurki-Suonio[[l986l:IIgnatius et al.L 
1994). Here the equation of state is that of a relativistic fluid, 
p = pc 2 /3, where c is the speed of light. Thus, there is no baro- 
clinic term and no obvious source of vorticity. However, the rel- 
ativistic equation of state may be modified at small length scales, 
but it is not clear that this can facilitate significant vorticity pro- 
duction. 



6. Conclusions 

The present work has demonstrated that vorticity production is 
actually quite ubiquitous once there is rotation, shear, or baro- 
clinicity. This implies that the assumption of potential flows as 
a model for interstellar turbulence might be of academic inter- 
est and could only be realized under special conditions of weak 
forcing, weak rotation, and no shear. In galaxies, however, the 
shear and Coriolis number are well below unity, leaving only 
the baroclinic term as a viable candidate f or the production o f 
vorticity. This agrees with early work of Kor pi et al.l dl999b), 
who analyzed the production terms in sup ersonic, supernova - 
driven turbulence quantitatively; see also iGlasner et al.l (fl997h . 
who showed that on long enough time scales significant vortic- 
ity can also be produced for subsonic flows. We have also ob- 
served how vorticity is mainly produced close to shock front en- 
counters. This motivates a more detailed investigation of these 
zones as the next step in the study of vorticity generation in the 
interstellar medium. It should also be pointed out that the baro- 
clinic Jerm_corre^OTidsto the battery term in the induction equa - 
tion iKulsrud et al.1 (fl997): Brandenburg & Subramanian (2005). 
Thus, when studying the possibility of dynamo action, this bat- 
tery term provides an intrinsic and well defined seed for the dy- 
namo and should therefore be included as well. 
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